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^ We investigate the potential of TESLA and JLC/NLC electron-positron 

D ■ linear collider designs to observe diquarks produced resonantly in processes 

■ involving hard photons. 
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I. INTRODUCTION 

Diquarks are amongst those particles whose features make them worth to be studied 
carefully. Like leptoquarks and dileptons, those particles are naturally predicted by many 



On ! interesting extensions of the standard model as well as are invoked in other models in order 
^ I to solve some of the most puzzling problems found in particle physics. These particles, 
O ' also known as biquarks and duarks, carry baryon number B = ±2/3 and lepton number 

^ ' L = 0, are colored and interact with quarks but not with leptons, and have integer spin. 

Scalar diquarks are predicted in many Grand Unified Theories (with or without super- 
Qh! symmetry) [1] and models of composite particles [2] , are invoked as an essential ingredient 
^ I for a model that evades the strong CP problem [3] and in [4] a class of models containing 
scalar diquarks is considered that can generate large CP conserving and violating contri- 
butions to double Cabibbo suppressed decays without affecting — mixing. Vector 
diquarks are predicted by models of composite particles [2]. 

Diquark production was already analyzed for e"'"e~, ep, pp and pp [5-8] colliders. They 
are expected to be easily observable at LHC [8] through their resonant production followed 
by the decay to two jets. Scalar and vector diquarks with masses up to approximately 
10 TeV could be observed, depending on their couplings to the quarks. If relatively 
light diquarks are observed they can also be expected to be on the reach of the next 
generation linear colliders, whose center of mass energy ranges from 500 GeV up to 1 TeV 
in the most conservative designs, like TESLA [9] and JLC/NLC [10]. The best channel 
for their observation on these machines would be, in principle, the pair production of 
diquarks. However, this channel is only useful if the diquark masses are smaller than 
\/s/2. This is a great limitation per se, and becomes more problematic when the present 
experimental limits on Eq superstring inspired model scalar diquarks mass is taken into 
account, 290 < M^^ < 420 GeV [11]. 



*e-mail:gusso@fisica. ufpr.br 



In order to overcome this limitation wc study the possibihty of observing diquarks 
produced resonantly through the process depicted in Fig. 1. The resonant production of 
new particles through higher-order processes has already shown its usefulness as shown, 
for example, in [12] for the resonant production o dileptons. In this article we will focus 
on the production of scalar diquarks. 

This article is organized as follows. In Sect. II we introduce the effective Lagrangians 
relevant to our analysis, and other theoretical and computational tools required to de- 
termine the cross-section for the process depicted in Fig. 1. In Sect. Ill we analyze 
the relevant background. In Sect. IV the results for the significance of the signal over 
background are presented. Finally, we conclude in Sect. V. 

II. SIGNAL 

Each extension of the standard model containing scalar diquarks we referred to in the 
Introduction has its won Lagragangian coupling diquarks with quarks. Yet, other models 
can be created and contain other Lagrangians. To take into account both the existing 
models and the ones that can be devised in the future we consider the most general 
SU{3)c X SU{2)l X U{1)y invariant Lagrangian coupling scalar diquarks with the quarks 
of the standard model. The coupling of diquarks to other bosons and with itself are of 
no relevance here. Following Ref. [13], the possible scalar diquarks are SU{3) triplets and 
sextets (In what follows, i,j and k are 5*^7(3) indices while a and /3 are SU{2), and family 
indices are omitted) 

5*(3, 1,-1/3), 5'*(3, 1,-4/3), S"'{3,3,-l/3), 5'"*(3, 1, 2/3), 
5^^(6,1,1/3), 5'*^(6,3,l/3), S"'\6,l,4/3), S"''\6,l,-2/3), (1) 

where the SU{3)c x SU{2) lxU{1)y multiplicities are given in parenthesis. Their couplings 
to the quark singlets u and d, and doublets q are given by 

+XrS'^dRiul^ + X8S"^qL^,{Ter''q'kp, + \^S"'^UiuU% + X^oS""^dRdly 

(2) 

From (2) we can get the Feynman rules required to determine the cross-section for 
the process in Fig. 1. In a process like this, most of the particles in the final state are 
expected to be lost in the beam pipe, the two exceptions being the quarks resulting from 
diquark decay. Consequently, the relevant signal is a pt balanced coplanar pair of jets. 

The cross-section for the production of a pair of jets through the process in Fig. 1 is 
given to good accuracy by 

a{s)e+e-^2jets= [ dxf dy fpe{x) fp^ e{y)(^{s)qq' ^q"q'"5{xyS - s) (3) 

In this expression, fpe{x) corresponds to the probability density of finding a parton p 
carrying a fraction x of the energy of the electrons (or positrons) in the initial beams. 
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a(s)qqr^q"q"' coriesponds to the cross-section for the subprocess in which two on-shcll 
quarks q and q', with total energy v^I = ^/xys, interact resulting into a diquark that 
utterly decays to the quarks q" and q'" . The lower limits of integration Xmin and ymim 
guarantees the energy required for the production of quarks. Because we are going to 
consider the five light quarks as massless those limits are actually equal to zero. 

The probability fpe{x) results from the convolution of the spectrum of photons emitted 
by the initial electrons, f-ydx), with fjyy{x) the photonic parton distribution function 
(PDF) 



dz dyf^e{z)fp^{y)S{zy - x) 



1 dz 
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(4) 



The photon spectrum /^e results from the combination of two effects, namely, 
bremsstrahlung and beamstrahlung. The photon spectrum resulting from bremsstrahlung 
for an electron beam with energy is [14] 
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The photon spectrum that results from the electromagnetic interaction between the elec- 
tron and positron beams in the intersection region, is described by a more complicated 
expression [15] 
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The beamstrahlung parameter T that characterizes each linear collider is given by. 



T 



(7) 



with Ee = Te = acm/nT-e = 2.818 ■ 10^"*^^ the classical electron radius, a^, (Jy and 

the average size of the particle bunches, and N the number of particles per bunch. Other 
definitions in (6) are 
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and the average number of photons per electron irradiated during the electron and positron 
bunches interaction 
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(9) 
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Because the photon production by bremsstrahlung and bcamstrahhmg arc distinct phys- 
ical processes the final spectrum is given simply by fyei^) = fjT"^^{x) + /^^""''(x). In 
Fig. 2 we present f-ye{x) for TESLA operating at ^/s = 800 GeV. For both TESLA and 
JLC/NLC at ^/s — 500 GeV we obtain a similar spectrum, however for JLC/NLC at 1 
TeV, the beamstrahlung is significant up to x ~ 0.7 instead. 

The photonic PDFs were determined by several groups each using distinct physical 
models and different sets of experimental data [16]. We adopted for our analysis the 
GRV photonic PDFs [17], calculated at leading order. This choice is based on two facts: 
(i) according to [16] the GRV parameterization gives one of the best descriptions of the 
structure function F2 and (ii) this parameterization is valid up to a characteristic scale 
q2 _ iq6 GeV^. Considering the characteristic scale for the scattering processes equal 
to s/4 implies that it is valid up to ^/s ~ 2 TeV. Recently, another parameterization 
[18] was developed, based on an improved model and more data, which fits the data on 
F2 better than GRV parameterization. However, both parameterizations fits the data 
almost equally well and the parameterization in [18] presents the limitation of being valid 
only up to ~ 900 GeV, while we are interested in energies up to 1 TeV. It is also 
worth mentioning that for the relatively high values of (the characteristic scale at 
which the PDFs are evaluated) and x (the energy fraction carried by the parton) that 
will be important in our search for heavy diquarks, the various parameterizations predict 
quite similar parton distributions [18]. In Fig. 3 we show fpe{x) for quarks and gluons as 
expected for TESLA operating at ^/s ^ 800 GeV, for = (400GeV)2. Similar results 
hold for TESLA and JLC/NLC. 

The last ingredient present in (3) is the cross-section o'{s)qqr^q'iqi" . In order to ob- 
tain the correct contribution from the subprocess qq' — > q"q"' we have to consider the 
interference of the s-channel production of the diquark with the diagrams involving the 
gluon exchange in both t and u channels. Denoting these contributions as a, b and c. 
respectively, the differential cross-section for the signal generated by a scalar diquark of 
mass Ms coupling to quarks with strength A can be written as 

da da{a) ^ da{a x b) ^ da{a x c) 



dt dt dt dt 

fl\2^. ]\/r2\ f'l\2^. ^l-s) 



: ^ 1 



167rs[(s - M|)2 + M|r2] t[(s - M|)2 + M|r2] m[(s - M|)2 + M|r2] ' 

(10) 



The decay width is 



FXM^ (11) 
167r ^ ^ 



The factors fs^f's^fj and Fg depend upon which of the interactions presented in (2) is 
considered and incorporate statistical factors associated with the presence of identical 
fermions in the vertices [19]. The sum and average over colors are incorporated in those 
factors as well. F^ contains the sums over all possible final states for diquark decay. The 
cross-section (10) was derived assuming massless quarks, a fairly good assumption at the 
energies we are going to consider. 



4 



At this point we have the freedom to choose which of the interactions from Cg are 
worth to be analyzed. The idea is to choose a representative interaction whose results 
can be lately interpreted approximately in terms of the other interactions. For that aim 
it is helpful to observe in Fig. 3 that fp^ for u and c quarks are at least three times bigger 
than for the other quarks. Because the cross-section (3) involves the product of two /pe, 
we are led to conclude that processes involving u and c quarks are about ten times more 
likely to occur. For that reason we choose the term XgS"^^URiU'[j from Cs, which involves 
u and c quarks. We also consider Ag, that must be symmetric, diagonal in the flavor 
space. Consequently the initial and final states in the subprocess in which the diquark is 
produced involve two u or two c quarks. For comparison, if Ag where flavor democratic, 
we could have 12 instead of 4 different subprocesses. The factors fs, f's, f'J and Fs are 8, 
8/9, 8/9 and 4, respectively, for the interaction we have chosen, for both u and c quarks 
in the vertices. 



III. BACKGROUND 

The two-jet background receives contributions from the annihilation process e'^e" — > 
^, Z ^ qq and the hard two photon process e+e" —>■ e~^e~qq [20]. The two photon 
processes can be of type direct, once resolved and twice resolved. The last two processes, 
despite being higher order in a^, contribute significantly at lower energies. 

The cross-section for the background generated by the annihilation process is given 

by 

Cr{s)annih^Y. fee{x) fee{y)cr{s)e+ e--^gg S{xyS - s) , (12) 



where fee{x) is the electron (and positron) spectrum, the probability density that we get 
an electron with a fraction x of the initial electron energy in the beams, and the sum runs 
over all possible quarks in the final state. cr(s)e+e-^gg corresponds to the cross-section 
for the process e'^e~ — > 7, Z — > qq. fee results from the convolution of the beamstrahlung 
and the initial state radiation (ISR) spectra 



feeix) = [ ^-^fi!\z)fer (f ) ■ 



We use the following analytical expression for f\ 



ISR, 



[21], 



/if^(^) = f (1 - ^)'-' (1 + y^) - + 



where /3 is defined as. 



/3 = ^""^^ 
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(13) 



(14) 



(15) 



The spectrum f^l"-"" was evaluated numerically using CIRCE [22] a FORTRAN routine 
containing parameterizations for electron, positron and photon spectra resulting from 
beamstrahlung at the future e+e" colliders. 
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The cross-section for the background generated by two hard photon processes is cal- 
culated from 

'^is)2^ = [ [ dy fieix)fjeiy)cr{s)ij^2jetsS{xyS - S), (16) 

where the indexes i, j = 7 or p. The sum over i and j corresponds to the sum over direct, 
once and twice resolved processes and the sum over p remind us that we have to sum over 
subprocesses involving different initial and final partons. In total we have to consider one 
direct process, and 2 once resolved, and 8 twice resolved subprocesses. The cross-sections 
for all such processes can be found in [23]. The five lighter quarks are assumed massless 
throughout our analysis and the top quark, of course, is not considered as a source of 
background. 



IV. RESULTS 

In the previous sections we presented all the theoretical and computational tools re- 
quired for our analysis. We used them to determine the expected significance of the signal 
over background for collider parameters as given in [9], for TESLA operating at 500 GeV 
and 800 GeV, and as given in [24] for NLC operating at 500 GeV and 1 TeV. Because 
the designs for JLC and NLC tends to converge, specially with respect to the parameters 
relevant for the beamstrahlung effect, we can in fact assume that our results are valid 
for JLC/NLC design. We show in Table I the relevant parameters for the calculation of 
beamstrahlung. 

The significance of signal over background is defined as. 



where S is the number of expected two-jet events attributed to the resonant production of 
diquarks, and B is the number of two-jet events resulting from the standard model alone. 
The signal generated by a diquark of mass Ms and decay rate F is calculated integrating 
the cross-section in the two-jet invariant mass (M) interval Ms — T < M < Ms + F, 
which embraces approximately 95 % of the events around the resonance. For a realistic 
analysis of the contribution of the background we considered the energy resolution of the 
hadronic calorimeter. For TESLA detector it reads. 

The same energy resolution was assumed for JLC/NLC. The corresponding two-jet invari- 
ant mass resolution is given approximately (we assume it to be exact in our calculations) 
by 5M = 0.5\/M + 0.02828Af . The background B is calculated by integrating the cross- 
section in the range Ms - AM < M < Ms + AM, with AM = max(F, 6M). 

The significance (17), scales with the collider luminosity as \/Z. Consequently, high 
luminosities are desirable for a good signal over background discrimination. Presently, 
TESLA is expected to deliver the highest luminosities, 500 fb~^/year, while the luminosity 
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for JLC/NLC is planned to be in the range 100-200 fb~^/year. With those luminosities 
in mind we determined the significance for the following total integrated luminosities: 
L — 100, 300 and 500 fb~^ for both collider designs and additionally L = 2000 fb~^ for 
TESLA. 

In determining 175 we considered some cuts on the pseudo-rapidity rj and the transverse 
momentum pt of the jets. We worked as out for I77I < 1 and 2, and \pt\ > 15 GeV, for 
both jets. This cut on pt is about the minimum required to guarantee that the assumption 
of massless quarks is valid for all quarks, including the bottom quark, because it implies 
an adequate minimum energy for the jets. It also guarantees the vahdity of the PDF 
parameterization and eliminates considerable amount of background. As we discuss in 
more detail later, higher Pt cuts did not improve as significantly, while reducing the 
absolute number of signal events to undesirably low values. Our analysis is also limited to 
the search for diquarks with masses above 100 GeV. Lighter Diquarks could in principle 
be pair produced by LEP, which reached more than 200 GeV as its center of mass energy. 
Here we are interested in what LEP could not reach. 

In Fig. 5 we present S/\fB as a function of the diquark mass, M5, for TESLA 
operating at 500 GeV. The solid and dotted lines are for |?7| < 1 and 2, respectively. 
The eight upper curves, both solid and dotted, correspond to a coupling with strength 
given by q; = A^/47r = 0.1, while the eight lower curves are for a = 0.01. For each 
a and cut on r] there are four curves corresponding to the integrated luminosities C — 
100,300,500 and 2000 fb~^. The curves for a = 0.01 where scaled by a factor 0.1 in 
order to clearly separate them from the curves for a = 0.1. The straight dashed lines 
correspond to S/^/B = 3, the minimum requirement for signal observability. The choice 
a = 0.1 corresponds approximately to the maximum acceptable value to keep higher order 
contributions under control. The coupling a — 0.01 is approximately the minimum value 
for which the signal significance is worth to be analyzed. For smaller values of a the signal 
significance is much smaller than 3 for masses greater than \/s/2. Fig. 6 shows the same 
curves for TESLA operating at 800 GeV. The results for JLC/NLC operating at 500 GeV 
and 1 TeV are presented in Figs. 7 and 8 for £ = 100, 300 and 500 fb~^. 

V. FINAL DISCUSSION AND CONCLUSION 

As already mentioned at the Introduction the Tevatron collider has set limits on 
the masses of a class of scalar diquarks, 290 < M^^^ < 420 GeV which are expected 
to be approximately valid for other scalar diquarks. There also exist indirect bounds 
imposed on Yukawa coupling constants for quarks and diquarks. Stringent indirect bounds 
for couplings involving the second and third generations where set from the limits on 
branching ratio for the inclusive decay B — > Xg^ [25]. Looser bounds on couplings 
involving other generations come from the analysis of precision LEP data [26], neutron- 
antineutron oscillation [27] and rare nucleon and meson decay [28]. The looser bounds 
generally permit couplings of the order one as required for a = 0.1 (A ~ 1.25) the higher 
value we admitted in our analysis. 

It is exactly for relatively high values of the Yukawa coupling constant that we may 
expect the resonant production of diquarks to be useful. This conclusion can be draught 
from the curves for a = 0.1 in Figs. 5, 6, 7 and 8. For couplings smaller by just one order 
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of magnitude (a = 0.01) the resonant production does not produce any interesting result. 
Is is worth mentioning at this point that for a = 0.1 the number of signal events is always 
of order of a few tens or more for 175 > 3, resulting in good statistics. Another conclusion 
that can be draught is that the best signal over background relation in the parameter 
region satisfying (75 > 3 and Mg > •\/s/2 is obtained considering the cut \r]\ < 1. However 
we can note from the results for the signal significance that the curves for jryj < 1 arc not 
always above the corresponding curves for \ri\ < 2. This behavior for cTs can be attributed 
to the distinct angular distributions for the background generated through annihilation 
and for the background generated by two hard photon processes. The two hard photon 
processes prevail at energies below -\/s/2 while the annihilation prevails at higher energies. 
This partially explains the behavior of the curves. An additional aspect of the background 
has to be taken into account, the fact that once and twice resolved processes initiated 
by gluons, that alter significantly the angular distribution, are more important at lower 
energies. Altogether, we have the elements that imply the relatively complex behavior of 
the curves for signal significance. 

A final comment on the pt cuts. We have tested cuts on pt higher than 15 GeV, 
however they did not improve the signal over background relation significantly. This can 
be explained by the fact that the dependence on px for the signal is almost the same as 
the one for the background generated by hard two photon processes and the background 
from annihilation is almost constant with respect to pt- 

We conclude, considering the present analysis and the previous works on the search for 
diquarks at the next generation of linear colliders, that diquarks with masses above ^/s/2 
are unlikely to be seen at such machines, unless new production channels are devised and 
prove to be sufficiently sensitive to the presence of those particles. The only exception 
are diquarks with relatively strong Yukawa coupling constants, that means A ~ 1, with 
the additional requirement of coupling to u and c quarks. 
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FIG. 4. The probability density fee for TESLA operating at ^ = 800 GeV. 
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FIG. 5. Signal significance for TESLA operating at s/s = 500 GeV as a function of diquark 
mass, Ms- Solid and dotted lines correspond to \r]\ < 1 and 2, respectively. For all curves 
\pt\ > 15 GeV. 




12 



lo- 



rn 
> 



00 



1 



: .V 1 ' 1 ' 1 ' 1 ' 1 




1 1 1 1 1 . 

JLC/NLC : 
Vs = 500 GeV 









\%i , 1 



100 150 200 250 



300 350 400 



450 500 



FIG. 7. Signal significance for JLC/NLC operating at ^/s = 500 GeV as a function of 
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TABLES 



TABLE I. Collider parameters relevant for the calculation of beamstrahlung. N is the 
number of particles per bunch, crx,y,z arc the mean bunch dimensions (z is the beam axis), T is 
the beamstrahlung parameter and the average number of photons per electron. 





Collider 


TESLA 


TESLA 


JLC/NLC 


JLC/NLC 


Parameter 




500 GeV 


800 GeV 


500 GeV 


1 TeV 


A[10iO] 




2.0 


1.4 


0.75 


0.75 


<^x [nm] 




553 


391 


245 


190 


ay [nm] 




5 


2.8 


2.7 


2.1 






300 


300 


110 


110 


T 




0.053 


0.084 


0.122 


0.315 






1.44 


1.39 


1.16 


1.38 
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